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Abstract 

Baldwin-Lomax and k-e turbulence models were modified for use in Navier-Stokes numerical 
computations of Mach 2.9 supersonic turbulent-boundary-layer flows along compression ramps. 
The computational results of Reynolds shear stress profiles were compared with experimental 
data. The Baldwin-Lomax model was modified to account for the Reynolds shear stress 
amplification within the flow field. A hybrid k-e model with viscous sublayer turbulence 
treatment was constructed to predict the Reynolds shear stress profiles within the entire flow 
field. These modified turbulence models were effective for the computations of the surface 
pressure and the skin friction factor variations along an 8° ramp surface. The hybrid k -e 
model could improve the predictions of the Reynolds shear stress profile and the skin friction 
factor near the corner of a 16° ramp. 


1. INTRODUCTION 

Shock-wave/turbulent-boundary-layer interaction is an important phenomenon which occurs 
in high-speed flow fields. Delery, Marvin, and Reshotko (ref. 1) have reviewed the literature on 
shock-wave/boundary-layer interaction in transonic and supersonic speeds. In their review, they 
describe flow phenomena based on experimental observations and theoretical concepts. 
Numerical solutions of the compressible Navier-Stokes equations and energy equation have been 
used by many researchers to predict the surface pressure and skin friction of flow fields with 
shock-wave/turbulent-boundary-layer interaction. The Reynolds analogy is usually used to 
relate the turbulent thermal flux to the turbulent momentum flux. Thus, accurate modeling of 
the turbulent shear stress is essential to the success of the Navier-Stokes computations. 

This paper describes some of the turbulence modeling techniques used in the Navier-Stokes 
computations of two-dimensional supersonic flow with shock-wave/turbulent-boundary-layer 
interaction. Flow fields around 8° and 16° compression ramps imbedded in fully developed 
turbulent-boundary-layer flows were considered in this study. (See fig. 1.) The free-stream 
Mach number for the flow fields considered was 2.9. A time-dependent explicit Navier-Stokes 
computational code (ref. 2) was used for the present computations. The Baldwin-Lomax 
turbulence model (ref. 3) was used first to model the turbulent eddy viscosity. This turbulence 
model was modified so that prediction of the Reynolds shear stress measurements could be 
made. A compressible k-e turbulence model (ref. 4) was then used to model the turbulent eddy 
viscosity of the supersonic, flat-plate turbulent boundary layer. The high-Reynolds-number form 
of the k-e turbulence model and a viscous sublayer turbulence model (ref. 5) were combined 
to formulate a hybrid turbulence model for the flow along a compression ramp. The computed 
Reynolds shear stress and mean flow properties within the compression ramp flows are described 
and compared with experimental measurements (ref. 6). 
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skin friction factor, r w /0.5 p^ U w 
empirical constants 
Reynolds number, p 
mass-averaged velocity components 
frictional velocity, (r w / p w )°* 5 
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mass-averaged turbulent kinetic energy 
surface pressure 
boundary layer thickness 
Kronecker delta 

mass-averaged turbulent dissipation rate 
turbulent eddy viscosity 
dynamic viscosity 
mass-averaged density 

Reynolds shear stress, dY + dV / 5X) 

surface shear stress, /i w ( dU/ dY)|^ 
vorticity, (dU/dY) — ( dV/ dX) 


Subscripts: 

c 

s 

V 

w 

0 

00 


condition at ramp corner 
sublayer condition 
sublayer edge condition 
wall surface condition 
condition at X = 0 
free-stream condition 


3. TURBULENCE MODELS 

This section describes the Baldwin-Lomax and k-e models, which have been widely used 
to model the turbulence in the computations of various types of turbulent flows. 

3.1. Baldwin-Lomax Turbulence Model 

The Baldwin-Lomax turbulence model is a two-layer eddy viscosity model (ref. 3). The 
eddy viscosity of a two-dimensional boundary layer can be written as follows: 

In the inner layer, the eddy viscosity /J t>i is given by 

P t ,i = p [cjY (l - e- Y+ / 26 ) ]%! (1 


In the outer layer, the eddy viscosity fi t o is given by 
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F max is the maximum value of 
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and Y max is the value of Y at which 
occurs. The function F kleb is the Klebanoff intermittency factor given by 
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The quantity is the difference between the maximum and minimum velocities. 
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The \VU 2 + V 2 j min term in U ( y f was set to be zero in this study. 

The empirical constants C-, C Q , and C kleb appearing in equations (1) to (3) were previously 
determined (ref. 3) for constant-pressure boundary layers at transonic speeds. In this study, 
these empirical constants were adjusted to achieve better prediction of the Reynolds shear stress 
in the supersonic compression ramp flow field. Different sets of these empirical constants were 
found to model the turbulence eddy viscosity upstream and downstream of the compression 
corner. These new sets of empirical constants are (l) Cj = 0.41, C Q = 1 and C kJeb = 0.8 for the 
upstream modeling, X < X c , and (2) Cj = 0.45, C Q = 2, and C kleb = 0.8 for the downstream 
modeling, X > X c . 


3.2. Hybrid k-e turbulence model 

Nichols (ref. 4) extended an incompressible k-e turbulence model (ref. 7) to include the effect 
of compressibility. Based on the existing turbulence modeling technique (ref. 4), the following 
compressible k-e turbulence model was used to predict the turbulent eddy viscosity within the 
supersonic, flat-plate turbulent-boundary-layer flow. 

The turbulent eddy viscosity was defined as 

,, = ,C„(kV,)(l - (4) 

with = 0.09. 


The turbulent kinetic energy and the dissipation rate were described by using the following 
transport equations: 

£ ( P k) + _£ (puk) + _£ (pvk) =£. 
at ax dY ax 



+ ± 

(m + M t )— 

1 axj 

dY 

ayj 


+ P - pe 


2/xk 

Y 2 


( 5 ) 


3 



- (pe) + — ipVe) + JL (pVe) =. 

at ax aY 


d 

(m + — ) 

de 

+ 1 

f . \ de 

(M + — ) — 

ax 

1.3 

dX 

dY 

1.3 dYJ 


+ 1.35 1 P - f {l.Sfps * ** 
k k 


-0.5pU r Y/p 


-) 


( 6 ) 


where 
f = 1 
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Some additional terms (ref. 4) to account for the compressibility effect on the transport of the 
turbulent kinetic energy were not retained in equation (5). The production term P was defined 
as 


p = W 


au; 


+ 


aUj 


ax, ax. 


2 auk 


3 ax v 


- - P M y f 5 

3 1J | aXj 


Equations (4) to (7) were also solved for the turbulent eddy viscosity in the computations of 
the flow along the compression ramp. The source terms were found to be numerically stiff at 
locations near the compression corner. Therefore, a viscous sublayer model (ref. 5) was used to 
eliminate the stiffness problem. This viscous sublayer model and the high-Reynolds-number 
forms of the k-e turbulence model corresponding to equations (4) to (6) were combined to 
formulate a hybrid k-e turbulence model for the flow over the compression ramp. 

To implement the sublayer model, it was assumed that 

(1) Adjacent to the wall, the viscous sublayer thickness Y y is defined by 

( 8 ) 


v 

where 


20 ((../(>„) <“/U, 

C is an empirical constant discussed later. 


The turbulent kinetic energy, the 


dissipation rate, and the eddy viscosity at the sublayer edge are given by 

K = (uX' 5 ) 
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(9) 
( 10 ) 

/V = <Vv k v/*v 

(2) Within the viscous sublayer, 0 < Y < Y y , the turbulent kinetic energy and the dissipation 


rate are given by 


k, = Mp.) k,(Y,/Y v ) ! 

(12) 

and 


*8 = 2 (^w/^w)( k v/ Y y) 

(13) 


By substituting equations (9) and (10) into equation (11), the following eddy viscosity 
expression is established: 
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It was further assumed that equation (11) holds near the sublayer edge. By substituting 
equations (12) and (13) in equation (11) and taking the limit as Y approaches 
following eddy viscosity expression is found: 
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Therefore, the eddy viscosity satisfying equations (14) and (15) and the limiting condition at 
Y y has the following form: 

Ct,» = C.(Y)(p/Y*)(k>,) (16) 

C g (Y) is assumed to be linear in the direction normal to the surface, that is, 

c,(Y) = C„ + (c, v - c„)y/y, (11) 

By substituting C g (Y) into equation (16) and evaluating the eddy viscosity at the sublayer 
edge, the expression 

C., = 0 - 8 (<V c °' 5 ) (/>v//O 0 ' s < 18 > 

can be obtained. Thus, the expression for the eddy viscosity within the viscous sublayer 
becomes 


ft,. = {c„ + McX' 5 ) - cy y/y v } (,/ y *)k ]/., (is) 

The viscous sublayer thickness is an important parameter of this sublayer turbulence model. 
For a fully developed flat-plate turbulent-boundary-layer flow with Mach 3 free stream, 
experiments (ref. 8) have established the correlation 

Y y /6 = 7.92xl0 2 Re'°' 91 (20) 

among the boundary layer thickness, the Reynolds number, and the viscous sublayer thickness. 
Based on the present flat-plate turbulent-boundary-layer conditions (Fig. 1), this correlation 
gave a viscous sublayer thickness of 4.8x10 2 mm. With the same boundary-layer flow 
conditions, equation (8) predicted the sublayer thickness when C y has a value of approximately 
0.4. Therefore, a C y of 0.4 was assumed for the viscous-sublayer eddy viscosity expression 
given by equation (19). 

In this study, equations (12), (13), and (19) were used to calculate the turbulent kinetic 
energy, the dissipation rate, and the eddy viscosity within the attached viscous sublayer. The 
high-Reynolds-number forms of equations (4) to (6) were solved numerically for similar 
turbulence quantities at locations outside the viscous sublayer. 


4. COMPUTATIONAL METHOD 

For this study, a time-dependent explicit numerical code (ref. 2) was used. Initially, this code 
used only the Baldwin-Lomax eddy viscosity model. The code was modified by adding the 
hybrid k-e turbulence model to the numerical procedures. The conservation equations, the 
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coordinate transformations between the physical and the computational domains, and the Finite 
difference scheme used are described in reference 9. A computational domain of 76 in the X- 
direction and 26 in the Y-direction was used for the 8° compression ramp computation, and a 
domain of 7.5£ in the X-direction and 2.56 in the Y-direction were used for the 16° 
compression ramp computation. These domains were chosen so that the shock waves crossed 
the downstream outflow boundary. Typically, a 221x149 mesh size was used for the present 
computations. The grid configuration was packed near the surface and around the compression 
corner. There were 12 grid points within the viscous sublayer. 

First, the turbulent flat-plate boundary-layer flow properties were calculated. The l/7th- 
power-law velocity profile was used as the initial mean flow velocity. The initial temperature 
profile was determined from the velocity and temperature relation of a compressible laminar 
boundary layer along a heated plate (ref. 10). Nonslip boundary conditions and constant wall 
temperature were used for the surface boundary conditions. The incompressible flat-plate 
boundary-layer turbulent kinetic energy profile (ref. 11) was used as the initial condition. The 
corresponding dissipation rate was calculated by using s = pC^k 2 //i t , where /i t was obtained 
from the Baldwin-Lomax model. The flat-plate boundary-layer flow properties were then used 
as the initial conditions to start the compression ramp flow computations. The free-stream 
properties (ref. 6) were used as the far-field boundary conditions. 


5, RESULTS AND DISCUSSION 

In this section, the computational results for the flow properties, especially the Reynolds 
stresses, are described and compared with measurements of reference 6. 

5.1. Flat-plate turbulent-boundary-layer flow 

The computed Reynolds shear stress profiles with the standard and revised Baldwin-Lomax 
turbulence models and the low-Reynolds-number k-e two-equation model are plotted in Fig- 
ure 2. The standard Baldwin-Lomax model calculated larger Reynolds shear stress at locations 
away from the wall surface. By changing the value of Cy e ^ from 0.3 to 0.8, and thus changing 
the Klebanoff intermittency factor (eq. (3)), the Baldwin-Lomax model can predict the 
experimentally measured Reynolds shear stress levels. The low-Reynolds-number k-e 
turbulence model (eqs. (4) to (7)) predicted the Reynolds shear stress measurements within the 
entire boundary layer. Particularly noteworthy, this turbulence modeling technique could 
describe the near-wall Reynolds shear stress variation as shown by the experimental results. 
The present computations also predicted a skin friction factor of 0.00095, which agrees 

with the measurement (ref. 6). 

5.2. Eight-degree compression ramp flow 

The computed Reynolds shear stress distributions along the Y-direction at five different X/6 
locations of the 8° compression ramp are shown in Figure 3. The Baldwin-Lomax model 
predicted lower downstream shear stress than the measured values. The modified Baldwin- 
Lomax model calculated larger downstream Reynolds shear stress levels, which agrees with the 
measurement. The present computations predicted the Reynolds shear stress variation within 
the entire flow field when the hybrid k-e model with = 0.075 was used for turbulence 
modeling. Preliminary study indicated that the requirement of this C/i value was due to the 
grid configuration used in the computation. 

The computed skin friction factor and surface pressure variations along the surface of the 
compression ramp are shown in Figure 4. The experimental surface pressure distribution can 
be predicted by using either the original or the modified Baldwin-Lomax model. The 
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computations indicate a higher level of the skin friction downstream of the compression corner 
when the modified Baldwin-Lomax model was used for turbulence modeling. This skin friction 
distribution agrees better with the measurements. The computation using the hybrid k-g 
turbulence model also predicted approximately the experimental surface pressure and skin 
friction distributions. 

5.3. Sixteen-degree compression ramp flow 

The computed Reynolds shear stress variations in the Y-direction at Five different X/6 
locations downstream of the corner are shown in Figure 5. Both the modified Baldwin-Lomax 
model and the hybrid k-g model could predict the Reynolds shear stress proFile only at the far 
downstream location (X/6 = 6.5). However, the hybrid k-g turbulence model prediction seems 
to show an improved prediction of the Reynolds stress proFiles near the compression corner. 

The computed skin friction factor and surface pressure variations along the 16° compression 
ramp surface are compared with the experimental measurements in Figure 6. The computations, 
using the Baldwin-Lomax model, closely predicted the surface pressure measurements. 
However, the computational results of the skin friction factor are in poor agreement with the 
experimental data. When the hybrid k-g model with C = 0.09 was used for turbulence 
modeling, the computation accurately predicted the downstream skin friction variation, and a 
lower level of the downstream surface pressure was computed. A large compression angle (16°) 
might cause the viscous sublayer thickness to be reduced downstream of the compression corner. 
In order to model this behavior, a smaller sublayer thickness (80 percent of its upstream value) 
at the locations downstream of the compression corner was used in this part of the 
computations. Without the reduction in the sublayer thickness, the computations would predict 
a lower level of the downstream skin friction. 


6. CONCLUSIONS 

An application of turbulence modeling techniques used in the Navier-Stokes computations of 
two-dimensional, turbulent-boundary-layer flow with shock wave interaction was investigated. 
The supersonic flow fields due to 8° and 16° compression ramps imbedded in fully developed 
turbulent-boundary-layer flows were considered in this study. The free-stream Mach number 
of the flows considered was 2.9. Based on the results of this study, the following conclusions can 
be made: 

1. The Klebanoff intermittency factor of the Baldwin-Lomax turbulence model must be 
modified for the model to predict the Reynolds shear stress profile within the flat-plate 
supersonic turbulent boundary layer. Some empirical constants of this turbulence model were 
also adjusted for the model to estimate the downstream Reynolds shear stress of a compression 
ramp flow. 

2. A hybrid k-g turbulence model, using viscous sublayer model for near-wall turbulence 
treatment, was constructed for the computation of supersonic boundary layer flow along 
compression ramps. The viscous sublayer thickness was an important parameter of the model. 
An experimental correlation for the sublayer thickness was used to resolve the effect of the 
sublayer thickness on the turbulence model. 

3. When the modiFied Baldwin-Lomax model or the hybrid k-g model was used in the 
computations of the flow over an 8° compression ramp, the computation could predict the 
experimental results of the Reynolds shear stress proFiles within the entire flow Field. The 
computational results of the surface pressure and the skin friction factor were found to agree 
well with their measured values. 
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4» The computations could predict only the Reynolds shear stress profile at far downstream 
locations when the modified Baldwin-Lomax model or the hybrid k-c model was used in the 
computations of the flow over 16° compression ramp. However, the hybrid k-e model 
improved the predictions of the Reynolds shear stress profiles and the surface skin friction 
factors immediately downstream of the compression corner. 
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Figure 2. — Boundary -layer-flow Reynolds shear stress 
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Figure 4. Surface pressure and skin friction factor of 0° compression ramp. 
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Figure 5. — Reynolds shear stress profiles of 16° compression ramp. 
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Figure 6. — Surface pressure and skin friction factor of 16° compression ramp. 
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